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We study delocalization transition in a one-dimensional Dirac fermion system with 
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random varying mass by using supersymmetric (SUSY) methods. In a previous paper, 
we calculated density of states and found that (quasi-)extended states near the band 
center are enhanced by nonlocal correlation of the random Dirac mass. Numerical 
studies support this conclusion. In this paper, we shall calculate localization length as 
a function of correlation length of the disorder. The result shows that the localization 
length is an increasing function of the correlation of the random mass. 
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1 Introduction 



Random disorderd system is one of the most interesting problem in condensed matter 
physics. Especially, localization phenomenon plays an important role in various prob- 
lems, e.g., qunatum Hall effect, transport properties of mesoscopic systems, quantum 
chaos, etc. In a previous paper (which will be refered to as paper I hereafter) |l[ , we 
studied random-mass Dirac fermion in one-spatial dimension, which is a low-energy 
effective field theory of random hopping tight-binding model, quantum spin chain, 
etc. This system has been studied extensively but in most of analytical studies only 
the 5-function type white noise limit is considered for technical reason fj, [|, |j. In 
most of realistic systems however, disorders have nonlocal correlations, and moreover 
sometimes correlation length of disorders is larger than typical length scale of the 
system. An important example is the quantum Hall state where magnetic length is 
smaller than averaged length scale of random potential. Another example is quan- 
tum spin chain with the ground state of dimer structure whose low-energy excitations 
are described by a Dirac fermion with random varying mass, i.e., the system under 
investigation. There impurities generate low-energy excitations. 

In paper I, we studied effects of nonlocal correlation of random mass, which is quite 
natural for application to realistic systems, and obtained density of states (DOS) as 
a function of correlation length of the disorder. The result indicates that (quasi- 
extended states near the band center are enhanced by nonlocal correlation of the 
random-mass variables. Recently we have performed numerical studies of the model 
and verified that analytical results obtained by supersymmetric (SUSY) methods are 
in good agreement with numerical calculations 0. 

In this paper, we shall study the same model and calculate localization length 
by using the SUSY methods. In Sect. 2, we shall briefly review the model and the 
SUSY methods. In Sect. 3, localization length is calculated by solving eigenvalue 
problem with respect to transfer 'Hamiltonian' for the spatial direction. Section 4 
is devoted for discussion. Physical meanings of the result of the localization length 
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and the DOS obtained in paper I will be discussed. Numerical studies also give us 
useful informations on the system. As an application of the results, we shall discuss 
low-energy properties of qunatum spin chains. 



2 Model and SUSY 

In this section, we shall review the model and SUSY methods mainly in order to fix 
our notations and make this paper self-contained. Reader who is familiar with the 
subjects can skip this section and go directly to Sect. 3. 

2.1 Model and Green's functions 

We shall study random system whose Hamiltonian is given by 

H c = - jdx ipyd^R - i\)\id x i\) L - im(x)(i/j\ l ip L - tp{ip R ) , 
= — Jdx ^hijj, 

h = -ia z d x + m(x)a y , V = (^l)*. (2.1) 
The random mass m(x) is decomposed into a uniform and random piece as 

m(x) = m + 0(x), (2.2) 

where [0] eras = and 

[<t>{x)<t>{y)]ens = 7^exp(-|x-y|/A), (2.3) 

with positive parameters g and A. The symbol [A] ens denotes the ensemble average 
of A over the randomness. It is easily verified 

ydx [(f)(x)(f)(y)] ens = g, 

and 

[(j)(x)(j)(y)] ens -> gS(x - y) as A -> 0. 



It is obvious that the parameter g controls magnitude of fluctuation and A is the 
correlation length of the disorder. 

Single fermion Green's function at energy u is defined as 



G af3 (x,y;iuj) = (x,a\ 



h — iuj 



\y,P), a,P = R,L, 



(2.4) 



where \x, a) is the normalized position eigenstate of fermion at x and chirality a. By 
functional integral, 

G a p(x,y;iu) = i^a^^y))^ 



— / V^Vij^ a (x)^l(y)e- s ^ 



dx ip\ih + uj)ip, 



S^p — 
Z\p = 

Ensemble averaged Green's function is obtained from (|2.4j ) as 
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(2.5) 



G a/3 (x-y;iuj) = (x,a 



h — iuj 

where the ensemble average is taken with respect to 4>(x) in Q2.1|) and (|2 
to ( |2.3| ). Introducing the bosonic superpartner £, 

G af3 (x - y; iu) = i(ip a (x)ip P (y)) • 



where 



with 



(Ah 



S= /dx 



The ensemble average can be converted into the functional integral form, 



(x)(f)(y)...] 



S r k 



V(f)((f)(x)(j)(y)...) exp(-S^), 
|dx l0(x)(-A 2 ^. 2 + l)0(x) 



(2.6) 
according 

(2.7) 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



From Eqs.( p.8| ), Q2.9| ) and ( |2.10| ), the expectation value of operator A is given by 

(A) s = J V^V^V£V^V(f) A exp(-(S + S+)). (2.12) 

The above total action is obviously invariant under SUSY transformation ip «-> £ and 
has such a form that the SUSY partners ip and £ couple to the 'dynamical' real scalar 
field cf) which is 'singlet' under SUSY transformation. 

Transfer Hamiltonian is obtained by regarding the spatial coordinate x as time, 
and then the system reduces to a quantum mechanical system with the two bosonic 
and one fermionic variables. By solving the Schrodinger equations in the quantum 
system, the Green's functions are calculated. 

2.2 Transfer Hamiltonian 

The transfer Hamiltonian is obtained by regarding the spatial coordinate x as time 
in the functional integral representation Q2.12 ). Then the system reduces to a quan- 



tum mechanical system. In Ref.||, the following canonical creation and annihilation 
operators are introduced corresponding to the functional integral variables, 

?P R ^F h ^l^F\, ^{^F b 

£ R ^B h &^B\, (2.13) 

Fermionic and bosonic spin operators are defined by, 

S = VaF, 

J= l -BoB, (2.14) 

where 

B = B ] a z . (2.15) 
The commutation relations followed by F a and B a is 

{F a ,Fl} = 6 a0 , (2.16) 
[B a ,B f j\=5 afl . (2.17) 



— # — * 

It is proved that S and J satisfy SU{2) and SU(1, 1) algebras, respectively, and they 
commute with SUSY charges Q and Q, 

Q = BF, Q = F j B, (2.18) 

which satisfy 

Q 2 = Q 2 = 0, (2.19) 

{Q,Q} = N = N B + N F , (2.20) 

with N B = BB and Np = F^F. Note that SUSY charges Q and Q are different 
from those in particle physics, where {Q, Q} = H with Hamiltonian H. F a and B a 
(o =|, |) part of the transfer Hamiltonian is written in terms of the spin operators. 
Transfer Hamiltonian of is also obtained from Eq. (|2.11| ). The system of <fi is nothing 



but a simple harmonic oscillator linearly coupled with the SUSY spin J = S + J. 
In terms of the spin operators and the canonical boson operators of the harmonic 
oscillator a, a) which correspond to 0, Hamiltonian of the system is given as, 

H = 2ooJ z + 2m J x + \%J x {a + a f ) + \{a ] a + \). (2.21) 

V A A 2 

Fermionic states of F a are specified by representations of SU{2). Similarly bosonic 
states of B a form multiplet of irreducible representations of SU(1, 1), which are spec- 
ified by total spin 

J 2 = {N B + 2N B )/4, N B = BB = B\B^ - B\B h (2.22) 

and ^-component of spin J z , i.e., 

J 2 \jn) = j(j + l)\jn), 



J z \jn) 



l + |2j + l| 



+ n 



\jn), (2.23) 



2 

where j = 0, ±1/2, ±1, ■ • ■ and n — 0, 1, • • -. 

Quantum states of the system are specified by the quantum numbers N B , N F , 
T = QQ, T = QQ, J z and N$ = a) a, which commute with each other. The operators 
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r and T satisfy identities, r 2 = T 2 = 0, T + T = N and IT = IT = 0. The state 
| N B , N F , T, T, J*, N^) is given by the direct product \N B , N F , T, T, J z ) <g> \m) H where 
N^\w)h = rn\m)H- In Ref. 0, the structure of \Nb, Np, T, f , J z ) is studied in detail. 
The state with Np = is 



|2j, 0, 0, N, (\N + 1| + l)/2 + n) = \jn) \vac), 



and the one with Np = 2 



\2j, 2, N, 0, (\N — 1| + l)/2 +n) = |jn)| IT), 



where the fermionic sector is given as 



t) = Fj\vac), 
l) = F}\vac), 
lV = FjF}\vac), 



(2.24) 

(2.25) 

(2.26) 
(2.27) 
(2.28) 



with | vac), the vacuum of the fermion F a . The states with Np = 1 can be constructed 
by acting Q and Q on the above states. For N ^ 0, we have 



\2j, 1, (N - \N\)/2, (N + \N\)/2, \N\/2 + n) 



n+\N\ 

I) 



n 



•b'n)|T) 5 



\\ 2n+ \N\ u '" n *' ' V 2n+ |JV| 
\2j, 1, (iV + |iV|)/2, (JV - |JV|)/2, |JV|/2 + n) 



n 



2n+ \N\ 



\jn)\i) + 



n+ \N\ 
\ 2n+\N\ 



\jn)\1). 



For = 0, however, these two sets of states coincide and are equal to 
I - M,0, 0,n) = - 1/2, n)\ |) + | - 1/2, n- 1)| T)), 



(2.29) 



(2.30) 



(2.31) 



where n > and | — 1/2,-1) = 0. Note that these are annihilated by both Q and 
Q. In order to complete the set of the quantum space, we introduce a set of states 
orthogonal to the states Eq. (|2.31|) , 



-l,l,*,*,n) = ^=(|-l/2,n)| I) 



l/2,n-l)| T», 



(2.32) 
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which are not eigenstates of T and T. 

In later discussions, we shall consider 'right' eigenstates of H, 

H\N B , N F , T, f , E) = E\N B , N F , T, f , E), (2.33) 

which can be expanded in a basis of appropriate eigenstates of J z and N^, 

\N B ,N F ,T,T,E)= ]T X N E B ' NF ' T ' f {J\mNB,N F ,T,T,J\N^). (2.34) 

Since Q and Q do not commute with N B , N F , T, f but do with H, they relate degen- 
erate eigenstates of H . For N ^ 0, there are two sets of SUSY doublet, 

Q\N,0,0,N,E) = f \N- 1,1,0, N,E), (2.35) 

Q\N -1,1,0, N,E) = -\N,0,0,N,E), (2.36) 

Jo 

Q\N,0,0,N,E) = Q\N -1,1,0, N, E) = 0, (2.37) 

and 

Q\N- 2,2,0, N,E) = ft\N - 1,1,0, N, E), (2.38) 

Q\N -1,1,0, N, E) = ^\N- 2,2,0, N,E), (2.39) 

Jo 

Q\N -2,2,0, N,E) = Q\N — 1, 1, 0, N, E) = 0, (2.40) 

where /o and /q are constants to be determined by normalization. For N = 0, the 
space spanned by | — 1, 1, *, *,n) is not closed under the operation of H. The best 
one can do in this sector is to find an eigenstate of H projected back onto the same 
sector, namely, 

if | - 1, 1, *, *, E) — E\ — 1, 1, *, *, E) + (2.41) 

with \ip) = J2 n i>n\ — 1) 1? 0, 0, n). Including the above state, one finds a superquadru- 
plet, 

Q\ - 1, 1, *, *, E) = /i|0, 0, 0, 0, E), (2.42) 
Q\-l,l,*,*,E) = f 2 \-2,2,0,0,E), (2.43) 



for iV ^ and 

Q|0,0,0,0,£)=/ 3 |-1,1,0,0,£), (2.44) 

Q\ - 2, 2, 0, 0,E) = f 4 \- 1, 1, 0, 0, E), (2.45) 
Q\ - 2, 2, 0, 0, E) = Q\0, 0, 0, 0, E) 

= Q\ - 1, 1, 0, 0, E) = Q\ - 1, 1, 0, 0, E) = 0, (2.46) 

for N = 0. Since the SUSY invariant 'vacuum' should be annihilated by Q and Q, 
the 'vacuum' has the eigenvalues Nb = —Np = 1 and r = f = 0, and therefore, it 
belongs to the subspace spanned by | — 1, 1, 0, 0, n) ® |m)#. 

We dedicate the rest of this section to review the 'left' eigenstate of H, SUSY 
invariant identity operators and supertraces 0. Please notice that the operator H 
is not hermitian, since J x is anti-hermitian. 'Left' eigenstate of H, therefore, is not 
identical with 'right' eigenstate. Since the operator 



transforms J x as 



U transforms 



U = e lnJ +1 * Nf (2.47) 

U ] J X U = - J x , (2.48) 

U ] HU = H\ (2.49) 

ipQU = (Q t , (2.50) 

rfQU = Q ] . (2.51) 

(2.52) 

Given a right eigenstate of H, \E) R , the state, 

\E) L = U^\E) R , (2.53) 

is left eigenstate of H, i.e., 

L (E\H = L (E\E. (2.54) 
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We normalize the states as 

L (E\E) R = 1. (2.55) 



This gives the constant f = f' Q = VN for A^ and h = f 3 , f 2 = / 4 for N = 0. 
For ^ 0, the SUSY invariant identity operators are 



l N ;N F =o,i = E I N, 0, 0, N, E) RL (N, 0, 0, N, E\ 

E ^ 

+\N — 1, 1, N, 0, E) RL (N - 1, 1, iV, 0,E\), (2.56) 



and 



1jv ; jv>=i,2 = 52(\N -1,1,0, N,E) R . L (N -1,1,0, N, E\ 

E V ~ 

+|JV - 2, 2, AT, 0, £> fi ■ l{N — 2, 2, A", 0, E^j. (2.57) 

Please notice [Q, 1] = [Q, 1] = 0.Q For A/" = 0, the normalization condition is more 
complicated since | — 1,1, 0, 0, E) R and | — 1,1,*,*, E) R are both orthogonal to their 
corresponding left eigenstates, i.e., 

L (-l, 1, 0, 0, E\ - 1, 1, 0, 0, E) R = 0, (2.58) 
L (-1,1,*,*,E\-1,1,*,*,E) R = 0. (2.59) 



We therefore require 



L {-1,1,0,0,E\-1,1,*,*,E) R =1, (2.60) 
L (-l,l, *,*,E\ -1,1, 0,0,E) R = 1, (2.61) 



with the usual normalization condition, 



L (0, 0,0,0, £10,0,0,0,^=1, (2.62) 
L (-2, 2, 0, 0, E\ - 2, 2, 0, 0, E) R = 1. (2.63) 



1 In order to have this property, the definition of U = e l7rjr ~ +I » Np is essential. This is the reason 
why we have chosen the different definition from that in Ref.|| 
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The corresponding identity operator for N = is therefore 

1at=o = E(|-2' 2 > °> °> E )R ■ l(~2, 2, 0, 0, E\ + |0, 0, 0, 0, E) R ■ L (0, 0, 0, 0, E\ 

E ^ 

+| -1,1, *,*,E) R - L (-l, 1,0,0, E\ + | -1,1,0, 0, J B) i? - L (-l,l, 

(2.64) 

In the following, we study the system in the long size limit, where the system size 
= L — > oo. In this case, the supertrace of an operator Oe~ LH is expected to reduce 
to an expectation value in the SUSY invariant vacuum, 

(O) = StrOe- LH - L (0|O|0) fl e- L ^, (2.65) 

where E b is the lowest eigenvalue of H. As explained above, the vacuum |0) R has 
eigenvalues N B = -N F = 1, F = f = 0. 



3 Localization length 

3.1 Eigenvalue problem 

The ensemble-averaged localization length is the parameter which measures how fast 
the avaraged one-particle fermion Green's function decays. The averaged fermion 
Green's function is given as follows, 

G aa {x-iu) = L (0\e xH F a e- xH Fl\0) R e-^. (3.1) 

Since the state Fl\0) R has the quantum number N B = —1, N F = 2, T = N = 1 and 

r = o, 

G aa {x-iu) = L(0\F a \ - 1, 2, 1, 0, E) R - L (-1, 2, 1, 0, E\Fl\0) R e~( E ~^ x e~^ 

E 

~ e't ) (3.2) 
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where £ w is 'imaginary-time' localization length and the state | — 1, 2, 1, 0, E)r satisfies 
the following equation, 

H\ - 1, 2, 1, 0, E) R = E\-1, 2, 1, 0, E) R . (3.3) 

It is obvious that the localization length is related with the eigenvalue E. We shall 
obtain the ginuine localization length £ e from £ w by the analytic continuation iuj — > e. 
In the following, we rescale the paremeters as uj — > go;, m — > gm and A — > ^. 

|— 1, 2, 1, 0, i£)_R is given by a linear conbination of the states \n)\ = | — 1, 2, 1, 0, n)^, 
and can be written as 



|-1,2, 1, 0, E) R = Xn,m\n)x\m) H , (3.4) 

n=0 

with coefficients Xn,m- Reminding that the wave function {0„ im } for |0)#, 

= E _ !' X ' °> °> n ) ® l m )^' 

has the form0, 

where 4>^ m (i = 1,2) are independent of A, and the 'energy eigenvalue' of \0)r is jr, 
we assume that Xn, m and are expanded as 

Xn,m — A 2 (xL,m ^Xn,m ^ Xn,m '"')■> (3-5) 

E = i- + + A^W + A 2 ^ 2 ! • ■ ■ , (3.6) 
2 A 

where Xn m an d -^^ are independent of A. Equation ( |3.3|) gives the Shrodinger equa- 
tion of Xn] 

£(#o + 4(^) 2 + ^)A*+f X jtWi 

= E^EA i+ *xa» 1 , (3.7) 
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where 

H = 2uJ z + 2m J x - 2(J X ) 2 . 

Comparing terms of each order of A leads us to the following series of equations, 
the 0(X~ l ) terms; 



E*S»i + \/J^ExL°L»i = 0, (3.8) 

n i n 



the 0(A°) terms; 



5>oxS»i 



the 0(A X ) terms; 



+4CT) 2 £ xg>>i + 2^+1 J x E tn\m + i\n)i 

n n 

+™EXninWl + 2 V / ^ExL 1 ,Ll^)l 
n n 

= ^ M Exg»i, (3.9) 

n 

E^oxJLl^i 

n 

+KJ x f E x£i»i + 2V^ttj* E x!!,U»i 

n n 

+^EXn,Ll^)l + 2 V / ^Ex5nl™)l 

n n 

= ^E4Un> 1 +^i;xg l »i, (3.10) 



and so on. It is easily seen that Eq. (|3~8|) and Eq. (|3.9| ) for m = determine Xh m - 
Equation( |3.9| ) for m = is identical with the equation of \ n which is considered by 
Balents and Fisher 0. In a similar way, Eq. (|3.9|) and Eq.( p.lO|) for m = determine 



X W and £?W. 



3.2 Solution in O(A ) 

The solution xf} m and E^ are given by Balents and Fisher ||. We review the 
derivation of them again for completeness. Equation to be solved is 

(2uJ z - AMuJ* + A{uJ*f)) ExE'oHi = £ M Exi 0l o|n)i, (3.11) 

n n 
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with M = 2rriQ. Using the representation of J z and J x for |n)i, 

J z \n) 1 = (n + ^)\n} 1 (3.12) 
J*\n) x = ^^|n + l) x - - l) l5 (3.13) 



2 1 ' 2 



(3.14) 



we have the descrete equation of Xn, , 

2^(n + IjxJJ, + 2M(n X i I 1)0 - (n + l)xl°li,o) - (n + 2)(n + l)xL°i 2 , 

+ (2n 2 + 2n + l)^ - n(n - l)xi°I 2 ,o = E^ Q . (3-15) 

Since multiplying a solution for M by (— l) n yields a solution for — M, we expect 

x!?o = c 3 X [( V, M, E) + c 4 (-ir X [0] (n, -M, E), (3.16) 

where x{ n i M, E) is obtained by solving the continuum equation derived from Eq.( |3.15| ), 

2un- 2M{2n— + 1) - (2n— + l) 2 ) X [0] { n , M , E ) = E l ° ] X l ° ] {n, M, E). (3.17) 



dn dn 
Changing variables n to z = Inn, Eq.( [3.17|) gives 



2— + l) 2 - 2M(2— + 1) + 2ue z )<5> E (z, M) = E [0] $ E (z, M), (3.18) 



dz dz 

where M) = x [0] ( n i M -, E )- 

Though it is possible to solve Eq. (|3.18| ) exactly, we use the hard- wall approxi- 
mation for simplicity, i.e., the 'potential' term ue z is neglected and the boundary 
condition lno;|, M) = is imposed instead. Then the solution is obtained as 

$^,M) = e-^sin^^, (3.19) 



with (3 = y/E — M 2 and z u = \ \nu\. In order to obtain a solution for xho> namely, 
to fix the constants C3 and C4 in Eq. (|3.16|) , additional conditions have to be imposed. 
Since neither ^t°l(n, M, E) = 1 nor (— l) n are, unfortunately, solutions in the limit 
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u) = M = 0, we are not able to use similar boundary conditions used in paper I in 
which we obtained the DOS. 

Our strategy is to compare the solution (|3.19|) with the one obtained by solving 
the discrete equation ( |3.15| ) by neglecting the term uj z . Since the other terms in the 
Hamiltonian except uj z depend on only J x , we consider eigenstate \a) of J x \ 

J x \a) = a\a). (3.20) 

This state is the eigenstate of H with u = 0, because 

H\a) = (4Ma - 4a 2 )\a) = E\a). 

Expanding \a) in terms of |n)i, 

\a) = Y,M<*)\n)u (3-21) 

n 

we find XnQ = C +V ; n( a +) + c_^ n (a_) with a± = M ^ . This approximation is valid 
for n 1/uj. Eq. (|3.21|) leads us to the following equation of ijj n (a), 

(n + l)ip n+ i(a) - n^n-i(a) = -2ai[> n (a). (3.22) 

Equation ( |3.22|) can be solved easily by introducing the function ip(w, a) = J2n " l Pn{<^)w n . 
By using the identity for a(w) = J2 n < 7 n'w n , 

1 Id 
2^(( n + 1 ) (J ™+i -na n ^)w n = -((l-w 2 )—-w)a(w) 

= D w a(w), (3.23) 



the difference equation ( |3.22|) is converted into the differential equation, 



d 

((1 - w 2 )- w + 2a)i>(w, a) = 0, 



dw 

and solution to this equation is obtained as 



ip(w, a) = c(l - w)-2 +a (l + w)' 
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with c being a constant. Thus we have 

c f dw (1 — w)~^ +a (l + w)~z 

Vn{Ot 



n\ Jc 2iri w n+1 

= fTTn E - 1 ) l T , {--l + n-a)T(- + l + a . 3.24 

r(n+l) 7r ^ 2 2 

The contour C is the circle whose center is located at the origin and radius is small 
enough. Dominant terms in Eq.( |3.24| ) are I = and/or / = n for large n. Using the 
fomula of the Gamma function, 

lim ^^±^ = 1, (3.25) 
n->oo Y(n)n a 

we find that, for large n, Eq.( |3.24[) can be approximated as, 

^(a) = n -^-a r (2+ a ) + (_ 1 )n w -|+a £(iZ^ > ( 3 . 2 6) 
'7T \/7T 



where we have fixed c = ^ without loss of generality. As we shall see, in the limit 

COS 7TO ° J ' 

M — > and u — > 0, a-t tend to vanish, and thus r(~ ± at) tend to y/n. Therefore we 
have 

X ® Q = e-^^c+e"^ 2 + c_e^) + (-l^e'^ic+e^ 2 + c.e"^ 2 ). (3.27) 



We shall compare Eq.( p.27|) with Eq. fl3.16Q . These two solutions should coincide 



with each other for 1/u -C n <C 1. This matching condition gives 

c + = c_ = C3/2 = C4/2, /3zu/tt = an even integer, 
c + = — c_ = c 3 /2z = — c 4 /2i, Pz^/tt = an odd integer (3.28) 

From the above condition on Pz^/n, E^ is determined. So much is the review for 
the zeroth order solutions. 

We summarize the solutions an d the eigenvalues E^ of Eq.( |3.11| ) as follows; 

xi ] o = c 3 (e-— z sin— (^-^) + (-ir + V— ^sin— (^-^)), (3.29) 
E ¥\ = m 2 + (—)\ (3.30) 



with an integer k > 1. The constant C3 is determined to be ^J 1 ^ 2 by the normalization 
condition. By solving Eq. fl3.8|), Xn.L( m > 0) i s determined in terms of Xnfi- 
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3.3 Solution in O(A) 

Now let us turn to , which gives the first-order correction of the localization length. 
From Eq. ( |3.9| ) and Eq. ( |3.10D for m = 0, we have, 

(Ho - E®) ExWjnh = --±-J*(H - £ [ °V*ExS>>i + £ x£l>>i. 

n III/ ~r X n n 

(3.31) 

Since E^ does not depend on m, it is sufficient to solve Eq. ( |3.31 ) for m = in order 
to obtain We employ a similar strategy used above, namely, we shall obtain 

solutions by both the hard-wall approximation to the continuum equation and the 
approximation which neglects the term uoJ z in H Q in the original descrete equation. 
The solutions depend on the constant E^K We obtain E^ by matching the solutions 
obtained by these methods in the region where the above two approximations are 
both legitimate. Using the representaion of J x and J z , Eqs. fl3.12"|) and ( |3.31| ) give 

(-(2^- + l) 2 - 2M(2— + 1) + 2ue z - E®)v E (z, M) 
v oz oz ' 

= -Aue z {2^- + l)$ E (z, M) + E^ E (z, M), (3.32) 
in the continuum approximation with n = e z and 

* E (z, M) = x ll] (e 2 , M) + 4(2^ + l) 2 <5> E (z, M), 

where 

x!!!o = c 3X [1] (n, M) + c 4 (-l) n X [1] (n, -M). 

In the hard-wall approximation, we neglect the term ue z on the left hand side of 

Eq.( |3.32| ), and we impose the boundary condition ^ E (z u! ,M) = instead. Making 

use of Eq. Q3~TP| ) and assuming ty E (z, M) = e~^~ z r](z, M), 

f) 2 8 M ft Pi 

(-4— - (3 2 )rj(z, M) = -8ue z (— - -) sin |(z - zj) + sin P -(z - zj). (3.33) 

It is not so difficult to solve this equation, 

, »^ 2 f z Wo xJ P, \ M . (3 em . p 

r)[z,M) = — — / dx[2ue ( — cos — (x — z w ) sm — (x — z^)) sin— (x — z u 

i\ , ) pj ZiU \ v 2 2 V ; 2 2 V ;; 4 2 V 
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x (^sin-(x - Z w ) cos -(2 - z a 



0, , . ( 

cos - (x - z u ) sin -(2 - z u 





+e z u 



u0{e z 







1 + 



1 



0, 



1+0 2 
0, 



l + (3 2 ~ 2 



0, 



2 sin-(z-^)) 



z , 1 Z 9 / 



sin -(« - ^)) + (e 2 - e z ") cos -(z - zj)) 



1 + 



(3 1/5 
— 4- ((2 - z«) cos - (2 - z w ) - sin -(z - zj)), 

apart from the general solution, sin \ [z — z u ), of the homogeneous equation. 



(3.34) 



On the other hand, the descrete equation ( [3.31D with ujJ z = is 

(4MJT - 4(^) 2 - E®) J2x l So\nh = £ X%\n)i- (3-35) 

n n 

In order to solve Eq. ( |3.35| ), we use Eqs.( |3.12"D and ( |3.23| ) as in Eq.( |3.24|) , and we have 
the following equation, 

(AMD W - AD 2 W - E®)x(w) = E^ X {w), (3.36) 

with x(w) = En Xn.W" ano - x( w ) = EnxlV 11, Decomposing each power of w in 
x(w), solution of Eq. (|3.36| ) gives us solution of Eq. (|3.35|) , Xn o' s - Equation (|3.36|) can 
be solved exactly. The solution is given by 

a + )ip(w, a-) - a + )ip(v, a_) / 



X (w) = 



dv 



(3.37) 



$(v, a + )ip'(v, a J) — V>'(f, a + )tp(v, a_) ^ (1 — t> 2 ) 2 ' 
This can be verified by directly substituting ( |3.37| ) into Eq.( |3.36j ). Performing the 
integral in ( |3.37|) and using the identity, 



$(v, a + )^}'{v , a J) — $(v, a + )ip(v, a_) = —2if3- ^(v, a + )ip(v, a_), (3.38) 



we have 



dv 



1-v 2 



$(v,ot + ) ${v,aJ) l-v< 



c+ I dv( ($(w, a+) - a+ % #(w, a_)) 



1 - 



+ 



),(3.39) 



where some factor in Eq.( |3.39|) can be simplified by 

ip(v, aq_) /l — v\ 4/3 



By using the identity, 



(3.40) 



1- V y0 l l /1-wyP 

dv\ = + const. 

1 + vJ 1-v 2 2i(3\l + wJ 

w , 1 1 . ,w + 1, . . 

c(f r = -In + const., (3-41) 

1 — it 2 w — 1 

we can calculate the above integral exactly and we have 

EK w - 1 

xM = 7777 In I — -~r\{c+^{w,a+) - c-ip(w, a_)), (3.42) 
4«p wj + 1 

up to the general solution x(u>). We hence obtain the exact form of Xn}o as 

V [H I dw xW 

Xn '° Jc 2m w n + l 
EK l 

= 4^E^7((- i r'- 1 )(^(« + )-c_^(«-)). (3.43) 

In order to compare the above result with the solution obtained by the hard-wall 
approximation, we approximate Xn]o as follows, 

w g[11 y 1 (( 1 1 . 

~ Ai(3 f^n-l^ C+ W{l,\-a + ) °~ IB{1, \ - a J) } 

- { - inC UB(,,l +a+ ) - C - ,B(,,l + aJ > < 3 ' 44 » 
where we have picked out the leading terms in ipi(a±), i.e., 

- JwA V + W -mn L v (3 ' 45) 

From the fact that — ^ \ , , has its maximum at I = n — 1 and 1 . , is a slowly 



varying function for large /, we find 



n— 1 j n— 1 j 



V— - ~ T 

f^n-l IB(1, § ± a) ^ n-l(n- l)B(n - 1, § ± a) 

/1 . T(n — |±q;) 
~ (Inn + 7) 7 r 2 - 

1 "r(n)r(f ±a) 

- (lnn + 7) ^! 1 , , n-5 ±a (3.46) 
1 (^ =■= a J 
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for large n, where 7 is Euler's gamma constant. We therefore have the first order 
correction of Xn,o as 



/ l + M _P §_ 1-M £ \ 

Xn,o - -T^(lnn + 7)(^n 2 (c + n 2 _ c _ n2 ) + (_l) n 2 ( c _ n 2 _ c+n2 )j 



E % + 7) (e^ 2 cos hz- zj) + (-l^+V 1 ^ 2 cos |(z - ^)) . 



In the limit <C 2; <C z w or c^e 2 <C 1, the solution in the hard- wall approximation 
is reduced to the form 

i+M E [1] 3 E [1] 3 

e-^( -jpz - ^(1 + M) — —z w ) cos P -{z - zj) + (-1)»(M - -M), (3.47) 

which should be compared with the solution obtained from the difference equation 
We hence obtain the 'energy eigenvalue' in 0(A), 

^--^■TW- {l + M) - {3A9) 

In the limit M — > and uo — > 0, we can neglect the higher-order terms in /3 and we 
have 

~ = _i_A; 2 . (3.50) 

Thus the energy eigenvalue is 

Since the one-particle Green's function has the form J2k f(k) exp(— E^x) with cer- 
tain calculable functions f(k), the 'imaginary-time' localization length £ w is the inverse 
of the coefficient of k 2 in E^, namely, 

r 2 1 r 2 A\ 
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The genuine localization length £ e is the real part of £ w after performing analytic 
continuation iu — > e + iu, and then Eq. (|3.52|) leads to 

lnfl 2 41 In ^| 

g \ TC Z 7T Z 

where we have recovered the constant q and use the fact that In - ~ In for small 

y 9 2 9 



lJlnfl 2 4|lnfk 
^-(M^ + ^^r^ , (3.53) 



e. Equation (|3.53|) shows that the localization length £ e is a increasing function of 
(gX). This means that the nonlocally-correlated disorder enhances the derealization. 
Moreover the number of states with energy below e is 

N(e) = fdep{e) 
Jo 

= -(- rr — — - ) + (higher-order terms), (3.54) 

2 ^ In 77- In 7T \ ; 

and this shows that the product of the localization length and the number of states 
is constant up to the first order of (gX), 

LN{e) = ^- 2 +0{{gXf). (3.55) 

Physical meaning of this result will be discussed in Sect. 4 from microscopic point of 
view (see also Ref. 0). 



4 Discussion 

In this paper, we have studied how the nonlocally correlated disorder affects the the 
mean localization length £ e in the random-mass Dirac fermions by making use of 
SUSY methods. (In fact this is a low-energy effective model of the random hopping 
tight binding model ||.) We found that £ e is an increasing function of the correlation 
length A. This result is consistent with the behavior of the DOS p(e) obtained in 
the previous paper QXJ] and the numerical studies in Ref. ||. There exists a simple 
intuitive picture of the result: For the case of larger correlation length A, average 
distance between kinks is longer and also average height of kinks is lower. Therefore, 
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finite correlation of disorders hinders the scattering of fermions by disorders and 
enhances the derealization of fermions. Though this result is obtained in the specific 
one-dimensional model, we expect that similar behaviours hold in general random- 
disordered systems. 

Let us consider relationship between the localization length £ e and the number of 
states N(e). Our SUSY calculation suggests that the following relation 

const. 

«■ = ~m (41) 

holds (up to 0(g\)). This can be understood as below. As is shown by the numerical 
calculations of wave functions of random-mass Dirac fermions ||], fermions in the 
random potential localize within an interval between adjacent nodes. According to 
'one-dimensional node counting theorem' [|j], the number of nodes A^(e) is equal to 
LN(e) — 1. Let us suppose that the position of each node Xiii = 1,2, ■ ■ ■ , No) is 
randomly distributed under the constraint < x\ < ■ • • < Xi < x i+ \ < ■ ■ ■ < x^ d < L. 
The expectation value (xj) is then given as 

1 f L dx . Na-l 



N k=1 Jo L l=1 



^=n / -f n o(x l+1 - Xl ). (4.3) 



with the normalization constant 

N * ^dxk N *- 1 
k=i uu L i=i 

Therefore it is found that the expectation value of the distance between adjacent 
nodes is (xj+i) — (x^ = 4. According to this estimation, we roughly understand that 
the localization length is to be inversely proportional to the number of states iV(e). 

In usual cases, it is not easy to obtain the localization length accurately by numer- 
ical calculations. However very recently, a very useful method for that was proposed, 
i.e., the non-Hermitian extension by introduction of an imaginary vector potential 0. 
Then we hope that the localization length in the present model is also calculated 
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numerically by that method, and the analytical expression, which we obtained in this 
paper, is compared with numerical calculations. This is under study and results will 
be reported in a future publication. 

While the results in this paper are obtained for the (effective) random hopping 
tight binding (RHTB) model, they give some important implications for closely re- 
lated model — spin-Peierls model @. Since the z-component of spin in the spin- 
Peierls model corresponds to the presence or absence of the tight-binding fermion, 
the Green's function in the RHTB model can be considered as a correlation function 
of spins in the spin-Peierls model. In fact, the single- fermion Green's function of the 
RHTB model, 

poo 

/ dte Kt (v\c n (t)cl(0)\v), 
Jo 

coincides with the spin-spin correlation function 

/ dte tet (v\S-(n,t)S + (n',0)\v), 
Jo 

where S (n, 0) are the spin- up and spin-down operators at site n, and S~(n,t) is 
given by 

S~(n,t) = e inspt S~{n,0)e- iHspt 

with the spin-Peierles Hamiltonian Tisp- The operators S + (n, 0) and c n are related 
with each other by the Jordan- Wigner transformation. In the spin systems, the mean 
localization length £ e , which is calculated in this paper, can be considered as the 
mean correlation length of the spins. We therefore find that the spin-spin correlation 
is enhanced by the suppression of randomness of disorders, because the suppression 
enhances the (quasi-) extended states near the band center e = 0. 

Let us compare the above result with the experimental observation in Ref.||. In 
Ref. ||, a phase diagram of Cui_ x Zna.Ge03 was obtained by the measurements of the 
magnetic susceptibility. The phase diagram shows that the antiferromagnetic (AF) 
order is enhanced by decrease of the impurity concentration at higher-doped region 
(x > 0.03). (Note that the spin-Peierls state collapses at around x = 0.030.) We 
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expect that the white-noise limit A = describes the system at around x ~ 0.5, i.e., 
extremely high-doped region. As increasing A, we approach the lower-doped region. 
Then as the calculations in this paper show, the extended low-energy excitations are 
enhanced as A — > large. This fact obviously means the enhancement of the AF order 
in the lower-doped region. This result is qualitatively in good agreement with the 
experiments. It is interesting to give more quantitative argument about this problem. 
This can be a useful check on our results. 
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